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ABSTRACT 

We present results of four-pointing Suzaku X-ray observations (total ^200 ks) of the intracluster 
medium (ICM) in the Abell 1835 galaxy cluster (kT ~ 8 keV, z = 0.253) out to the virial radius (r v i r 
~ 2.9 Mpc) and beyond. Faint X-ray emission from the ICM out to r v - lT is detected. The temperature 
gradually decreases with radius from ~8 keV in the inner region to ^2 keV at r v ; r . The entropy profile 
is shown to flatten beyond rsoo, in disagreement with the r dependence predicted from the accretion 
shock heating model. The thermal pressure profile in the range 0.3rsoo r ^ r v i r agrees well with 
that obtained from the stacked Sunyaev-Zcl'dovich effect observations with the Planck satellite. The 
hydrostatic mass profile in the cluster outskirts (?"5oo rS rS ?Vir 

) falls well short of the weak lensing one 
derived from Subaru/Suprime-Cam observations, showing an unphysical decrease with radius. The 
gas mass fraction at r V i r defined with the lensing total mass agrees with the cosmic baryon fraction 
from the WMAP 7-year data. All these results indicate, rather than the gas-clumping effect, that 
the bulk of the ICM in the cluster outskirts is far from hydrostatic equilibrium and infalling matter 
retained some of its kinetic energy. Finally, combining with our recent Suzaku and lensing analysis of 
Abell 1689, a cluster of similar mass, temperature, and redshift, we show that the cluster temperature 
distribution in the outskirts is significantly correlated with the galaxy density field in the surrounding 
large-scale environment at (l-2)r V j r . 

Subject headings: galaxies: clusters: individual (Abell 1835) — gravitational lensing: weak — inter- 
galactic medium — X-rays: galaxies: clusters 



1. INTRODUCTION 

Clusters of galaxies are the largest self-gravitating sys- 
tems in the universe, where thousands of galaxies and 
hot thin plasma (intracluster medium; ICM) are bound 
to the potential of the dark matter halo. Gravity of dark 
matter, which is the dominant mass component of clus- 
ters of galaxies, plays an important role in the structure 
formation and cluster evolution. According to the hier- 
archical structure formation scenario based on cold dark 
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matter (CDM) paradigm, less massive systems collapse 
first and then massive ones later. X-ray observables of 
the ICM properties keep original records of cluster evolu- 
tion. During the hierarchical formation, gas and galaxies 
in large-scale structure are falling on the clusters. Since 
cluster outskirts is located around the boundary of the 
cosmological environment, the gas in the outskirts would 
be significantly affected by structure formation. The 
cluster outskirts is, therefore, a good spot to refine the 
details of how the gas physics is involved in hierarchical 
clustering. It is, however, difficult to efficiently observe 
faint X-ray emission from cluster outskirts with Chandra 
and XMM-Newton because of their relatively high levels 
of instrumental background. 

Thanks to the low and stable particle background of 
the X-ray Imaging Spectrometer (XIS; Koyama et al. 
2007), Suzaku (Mitsuda et al. 2007) was able to unveil 
for the first time the ICM beyond rsoo , within which the 
mean cluster-mass density is 500 times the cosmic crit- 
ical density. Indeed, Suzaku's ability to probe the ICM 
out to the virial radius has been shown for a number 
of relaxed clusters (e.g. Fujita et al. 2008; George et al. 
2009; Reiprich et al. 2009; Bautz et al. 2009; Kawaharada 
et al. 2010; Hoshino et al. 2010; Simionescu et al. 2011; 
Akamatsu et al. 2011; Walker et al. 2012a,b; Sato et al. 
2012). One common feature is a flattening of the entropy 
profile beyond rsoo, contrary to the power- law prediction 
of the accretion shock heating model (Tozzi & Norman 
2001; Ponman et al. 2003; Voit et al. 2005). The entropy 
profiles, scaled with the average ICM temperature, are 
universal irrespective of cluster mass (Sato et al. 2012). 
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One possible explanation for the low entropy is devia- 
tions from hydrostatic equilibrium in the outskirts (e.g. 
Kawaharada ct al. 20f0; Sato et al. 20f 2). With simula- 
tions, Nagai & Lau (20f f) showed that beyond r2oo, gas 
clumping leads to an overestimation of the observed gas 
density. Simioncscu et al. (2011) interpreted, based on 
the results for Perseus cluster, that the entropy flatten- 
ing is a consequence of the gas density in the outskirts 
being overestimated due to gas-clumping. 

A gravitational lensing study is complementary to X- 
ray measurements, because lensing observables do not 
require any assumptions on the cluster dynamical states. 
Weak gravitational lensing analysis is a powerful tech- 
nique to measure the mass distribution from outside the 
core to the virial radius. The exquisite Subaru/Suprimc- 
Cam lensing data allows us to study properties of cluster 
mass distribution, thanks to its high image quality and 
wide ficld-of-view (e.g. Broadhurst et al. 2005; Okabe & 
Umetsu 2008; Umetsu & Broadhurst 2008; Okabe et al. 
2010a). Comparisons of X-ray observables with weak 
lensing mass allow us to conduct a powerful diagnostic 
of the ICM states, including a stringent test for hydro- 
static equilibrium (Okabe & Umetsu 2008; Kawaharada 
et al. 2010; Zhang et al. 2010; Okabe et al. 2010b). Kawa- 
harada et al. (2010) found in Abell 1689, incorporating 
Suzaku X-ray and lensing data, a large discrepancy be- 
tween hydrostatic equilibrium (H.E.) and lensing masses, 
especially discovered that H.E. mass significantly drops 
off in the outskirts (r > r5oo). 

Abell 1835 with an ICM temperature of —8 kcV is 
one of the luminous cool-core galaxy cluster. The X-ray 
properties of this cluster were measured within r$oo with 
XMM-Newton (Jia et al. 2004; Zhang et al. 2007) and 
with Chandra (Li et al. 2012). The temperature measure- 
ment in the outskirts were reported out to 9.'0 by XMM- 
Newton (Snowden et al. 2008) and to lO'.O for the western 
direction by Chandra (Bonamente et al. 2012), respec- 
tively. Okabe et al. (2010a) have conducted weak-lensing 
analysis of Subaru/Suprimc-Cam data to measure mass 
profile using the tangential distortion profile outside the 
core. Pereira et al. (2010) presented a complex velocity 
distribution, suggesting ongoing mass accretion associ- 
ated with smaller satellite systems and found that a third 
of Herschel sources are located in the southwest region. 
Morandi et al. (2012) presented full three-dimensional 
structure reconstructed from X-ray, Sunyaev-Zel'dovich 
(SZ) and strong lensing data available for the core region 
and discuss the non-thermal pressure with an extrapola- 
tion to r 2 oo- 

This paper reports the results of four Suzaku observa- 
tions of the Abell 1835 cluster out to the virial radius 
(r v ; r ~ 2.9 Mpc or 12.'0) and beyond. The Suzaku ob- 
servations and data reduction are described in Section 2. 
The spectral analysis to obtain radial profiles of temper- 
ature, electron density, and entropy is shown in Sec. 3. 
We discuss, in Sec. 4, a comparison of hydrostatic equi- 
librium (H.E.) and lensing masses, and gas mass frac- 
tion. We also compare thermal properties with those of 
other clusters, including stacked SZ pressure profile with 
Planck satellite. A statistical approach to investigate 
the correlation between temperature distribution in the 
outskirts and the large-scale structure derived from the 
SDSS photometric data, is conducted for a sample of two 
lensing clusters of Abell 1835 and Abell 1689, of which 



Suzaku data fully cover the whole region out to the virial 
radius. 

We use the Hubble constant Hq = 70 km s _1 Mpc -1 
(h = H I 100 km s _1 Mpc^ 1 = 0.7), assuming a flat 
universe with f2 m; o = 0.27 in this paper. The angular- 
size distance D A of the Abell 1835 is 819 Mpc h^ . This 
gives physical scale 1' = 237.9 kpc at the cluster red- 
shift z = 0.2532. We adopt the virial radius r v - n - = 2.89 
Mpc hjQ , within which the mean cluster- mass density 
is 112 times the cosmic critical density, determined by 
weak lensing analysis (Okabe et al. 2010a). The Galactic 
hydrogen column density nn of Abell 1835 is 2.04xl0 20 
cm~ 2 (Kalberla et al. 2005). The definition of solar abun- 
dance is taken from Loddcrs (2003), in which the solar Fc 
abundance relative to H is 2.95xl0 -5 . Errors are given 
at the 90% confidence level except as otherwise noted. 

2. OBSERVATIONS AND DATA REDUCTION 

2.1. Observations 

We performed four-pointing Suzaku observations of 
Abell 1835, named East, South, West, and North, in July 
2010 with exposure of ~50 ks for each pointing. The ob- 
servation log is summarized in Table 1. Figure 1 shows 
the XIS image of Abell 1835. The pointings were coordi- 
nated so that the X-ray emission centroid of Abell 1835 
was located at one corner of each pointing. The XIS mo- 
saic covered the ICM emission out to the virial radius 
(-2.9 Mpc or 12f0) and beyond. 

2.2. XIS Data Reduction 

We used only XIS data in this study. Three out of 
the four CCD chips were available in these observations: 
XIS0, XIS1, and XIS3. The XIS1 is a back-illuminated 
(BI) chip with high sensitivity in the soft X-ray energy 
range, while the XIS0 and XIS3 are front-illuuminated 
(FI). The instrument was operated in the normal clock- 
ing mode. We included the data formats of both 5x5 
and 3x3 editing modes in our analysis using xselect 
(Ver.2.4b). We used version 2.5.16.28 of the processed 
data screened with the standard filtering criteria 13 . In 
order not to reduce the exposure time, event screening 
with the cut-off rigidity (COR) was not performed in our 
data because we can estimate the non X-ray background 
(NXB) reasonably using data outside r v j r . The analy- 
sis was performed with HEAsoft ver 6.11 and CALDB 
2011-09-06. 

For each pointing (four azimuthal directions), we di- 
vided the field-of-view (FoV) to six concentric annu- 
lar regions centered on the X-ray emission centroid, 
(a,<5)=(14 h 01 m 01.865 s ,+02°52'35.48") in J2000 coordi- 
nates (Zhang et al. 2007), to obtain the temperature and 
electron density profiles. Inner and outer radii of the 
annular regions are 0^0-2^0, 2f0-4f0, 4f0-6.'0, 6.'0-9.'0, 
9.'0-12f0, and 12.'0-20f0 (Figure 1). The circular regions 
around 32 point sources were excluded from the analysis. 
Furthermore, we subtracted the contribution of luminous 
point sources outside the excluded regions. The details 
are described in Appendix A. In the O.'O— 2f0 region, all 
the spectra of the East, South, and North directions are 
unavailable because of calibration sources. 

13 http:/ /www. astro. isas.ac.jp/suzaku/process/v2changes/ 
critcria_xis . html 
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Fig. 1.— (left): NXB subtracted XIS image (0.5-8.0 keV) of Abell 1835. The XISO, XIS1, and XIS3 images were added on the sky 
coordinate, corrected for exposures, and smoothed by a 2-dimcnsional Gaussian with a = 16 pixels £S 17" (counts pixel -1 Ms -1 ). Here, the 
effect of vignetting was not corrected and regions where 55 Fe calibration sources are irradiated (Koyama et al. 2007) are excluded. Green 
circles indicate the regions used for spectrum analysis. Thick green circle shows the virial radius of Abell 1835 (r v j r ~ 2.9 Mpc or 12'.0). 
Small white circles are the excluded regions around point sources. Yellow boxes show the field-of-views (FoVs) of Suzaku Observations 
named East, South, West, and North, (right): XMM-Newton MOS1 + MOS2 image of Abell 1835 (0.5-2 keV). Background was not 
subtracted and vignetting was not corrected. 



TABLE 1 

Suzaku OBSERVATION LOG OF Abell 1835 



Observation (ID) Start a End a Exposure 

(ks) 



East (805037010) 2010/07/05 16:26:02 2010/07/07 02:55:18 49.4 

South (805038010) 2010/07/07 02:56:23 2010/07/08 11:12:11 45.6 

West (805039010) 2010/07/08 11:13:00 2010/07/09 23:39:12 53.7 

North (805040010) 2010/07/13 23:47:30 2010/07/15 09:36:15 48.8 



a Time is shown in UT. 



Redistribution matrix files (RMFs) of the XIS were 
produced by xisrmfgen version 2009-02-28. We 
generated two Ancillary response files (ARFs) by 
xissimarf gen version 2010-11-05 (Ishisaki et al. 2007), 
assuming uniform sky (circular region of 20' radius, 
here-after UNI-ARF) and surface brightness profile of 
Abell 1835, where used a /3- model image of 48f6 x 48f6 
with j3 = 0.55 and r c = 0fl92 based on the ROSATWRI 
result as the input X-ray image (Ota & Mitsuda 2004) 
(here-after (3-ARF). Using these ARFs, the normaliza- 
tion of the ICM component derived from the spectral 
fit for a given spatial region is that for the entire input 
region to calculate the ARFs. To derive the normaliza- 
tions of the ICM component for each spatial region, we 
multiplied the SOU RCE .RATIO. REG parameter from 
the xissimarf gen tool (see also e.g. Ishisaki et al. 2007; 
Sato et al. 2007) by the normalizations for the entire in- 
put region. We also included the effect of contaminations 
on the optical blocking filter (OBF) of the XISs in the 
ARFs. The NXB were estimated from the database of 
Suzaku night-earth observations using xisnxbgen version 
2010-08-22. 

3. SPECTRAL ANALYSIS AND RESULTS 



3.1. Spectral Fit 

We used the XSPEC vl2.7.0 package and AT0MDB v2.0.1 
for all spectral fitting. The NXB components were sub- 
tracted before the fit. To avoid systematic uncertainties 
in the background, we used energy ranges of 0.6-7.0 keV 
for the XISO, 0.5-5.0 keV for the XIS1, and 0.6-7.0 keV 
for the XIS3 in all the regions. In addition, we excluded 
energy band around the Si-K edge (1.82-1.84 keV), be- 
cause its response was not modeled correctly. We simul- 
taneously fitted all the spectra of the three detectors for 
the six annular regions toward four azimuthal directions 
by minimizing the total % 2 value. In this fit, relative nor- 
malizations between the three sensors were left free to 
compensate for the cross-calibration errors. The model 
for the spectral fit was an absorbed thin-thermal emission 
model represented by phabs x apec for the ICM emission 
of the cluster, added to the X-ray background (XRB) 
model. We employed the /3-ARF for the ICM compo- 
nent (see Section 2.2). The phabs component models the 
photoelectric absorption by the Milky Way, parameter- 
ized by the hydrogen column density that we fixed to the 
Galactic value of 2.04xl0 20 cm -2 (Kalberla et al. 2005). 
Even if we allowed the hydrogen column density for the 
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phabs model to vary or employ wabs model fixed at the 
Galactic column density in the direction of Abcll 1835 
in the spectral analysis, the resultant temperatures and 
electron densities of the ICM component are almost the 
same within ~3%. The apec is a thermal plasma model 
by Smith et al. (2001). For a given annulus, each pa- 
rameter of the ICM component for four azimuthal di- 
rections was assumed to have the same value. In the 
central regions, metal abundance of the ICM component 
was allowed to vary, while at r > 4.'0, we fixed the metal 
abundance of the ICM at 0.2. The redshift of the ICM 
component was fixed to 0.2532. 

In order to study the faint X-ray emission from the 
cluster outskirts, an accurate estimation of the XRB is 
vitally important. We fitted the spectra in the outermost 
annulus (12f0-20f0 region, which is outside r v ii) for the 
following three cases. Case- GAL: The XRB model in- 
cludes three components of the cosmic X-ray background 
(CXB), unabsorbed 0.1 keV Galactic emission (LHB; 
representing the local hot bubble and the solar wind 
charge exchange) and absorbed 0.3 keV Galactic emission 
(MWH1; representing the Milky Way halo; Yoshino et 
al. 2009). The normalization for the ICM flux is fixed to 
zero. Case-GAL+ICM: The XRB model was the same as 
the Case-GAL, but the temperature and normalization 
for the ICM component were left free. Case-GAL2: In 
addition to the Case-GAL, we added a absorbed 0.6 keV 
Galactic emission (MWH2; representing the Milky Way 
halo; Yoshino et al. 2009). This is because several blank 
fields observed with Suzaku contains the emission with 
0.6-0.8 keV (Yoshino et al. 2009). The normalization for 
the ICM flux was fixed to zero. In all cases, we assumed 
a power-law spectrum for the CXB with r = 1.4. In ad- 
dition, we modeled the LHB, MWH1, and MWH2 with 
apec model, where redshift and abundance were fixed at 
and unity, respectively. The temperatures of the LHB, 
MWH1, and MWH2 were fixed at 0.1 keV, 0.3 keV, and 
0.6 keV, respectively. We used the UNI-ARF for the 
XRB components, assuming that the XRB components 
have flat surface brightness (see Section 2.2). Normaliza- 
tions of the XRB components were also left free. As the 
XRB components, we adopted the model formula, phabs 
x ( powerlaw + apecyrwni ) + apecLHB for the Case- 
GAL and Case-GAL+ICM, and phabs x ( powerlaw + 
apecMwm + apecMWH2 ) + apecLHB for the Case-GAL2, 
respectively. Results of the spectral fit are shown in Sec- 
tion 3.2 and Section 3.4. 

3.2. Result of the XRB components 

Figure 2 shows the results of the spectral fit for the out- 
ermost annulus of 12f0— 20.'0 toward the East and West 
directions (opposite azimuthal directions). The best- 
fit parameters for XRB components and x 2 values are 
listed in Table 2. The x 2 values for the Case-GAL and 
Case-GAL2 are worse than the Case-GAL+ICM. We es- 
timate F-test probabilities for the Case-GAL and Case- 
GAL+ICM 'of ~1 x 10~ 7 and for the Case-GAL and Case- 
GAL2 of ~1 x 10~ 4 , respectively. The Case-GAL is there- 
fore not supported. The intensity of the Galactic emis- 
sions (LHB, MWH1) is somewhat higher than that of 
the typical Galactic emissions (Yoshino et al. 2009). The 
plausible cause of the higher intensity would be the fact 
that Abcll 1835 is located near the North Polar Spur. 
For the Case-GAL2, we refitted spectra with tempera- 



ture for 0.6 keV Galactic emission (MWH2) allowed to 
be a free parameter. We found that the \ 2 (2869) be- 
comes slightly better and the resultant MWH2 temper- 
ature increases to 0.92^09 keV. In the r < 12'.0 region, 
the best-fit ICM parameters were almost same as those 
from the Case-GAL+ICM. For example, the ICM tem- 
perature and electron density in the 9f0-12.'0 region in- 
creased by 8% and 10%, respectively, comparing to those 
for the Case-GAL+ICM. 

We investigated the validity of the CXB intensity ob- 
tained from the spectral fit. To estimate the amplitude 
of the CXB fluctuations, we scaled the fluctuations mea- 
sured from Ginga (Hayashida et al. 1989) to our flux limit 
and FoV area using the method of Hoshino et al. (2010). 
The fluctuation width is given by the following relation, 

^Suzaku ^Ginga f ^c, Suzaku j f ^c, Suzaku j ^\ 

ICXB ICXB V ^c, Ginga / \ S c , Ginga / 

where (<rs U zaku/^cxB) means the fractional CXB fluctu- 
ation width due to the statistical fluctuation of discrete 
source number in the FoV. Here, we adopt tJGinga/ Icxb = 
5%, with S c (Ginga: 6x 10~ 12 erg cm -2 s _1 ) representing 
the upper cut-off of the source flux, and O e (Ginga: 1.2 
deg 2 ) representing the effective solid angle of the detec- 
tor. We show the result, cr/Icxs, for each annular region 
in Table 3, where a is the standard deviation of the CXB 
intensity, Icxb ■ Table 3 shows f2 c (the solid angle of ob- 
served areas), Coverage (the coverage fraction of each 
annulus, which is the ratio of Q c to the total solid angle 
of the annulus), SOU RCE .RATIO. REG (the fraction 
of the simulated cluster photons that fall in the region 
compared with the total photons generated in the entire 
simulated cluster), and a /Icxb (the CXB fluctuation due 
to unresolved point sources). For all directions, a /Icxb 
values in the 9f0-12.'0 and 12.'0-20f0 regions are about 
6.1% and 3.3%, respectively. u/Icxb value for each di- 
rection is higher by a factor of ~2. These ranges are 
consistent with those obtained by the method of Bautz 
et al. (2009) using parameters derived by Moretti et al. 
(2003) studied ROSAT, Chandra, and XMM-Newton ob- 
servations. 

The best-fit parameter of CXB surface brightness (af- 
ter subtraction of point sources brighter than 2xl0~ 14 
erg cm" 2 s" 1 in 2-10 keV band) is 5.17±g;|| x 10~ 8 erg 
cm~ 2 s _1 sr _1 for the Case-GAL+ICM which agrees 
with Hoshino et al. (2010) of 4.73±££j xl0~ 8 erg cm~ 2 
s -1 sr _1 , within statistical errors taking into account the 
CXB fluctuation. Using the same threshold, Scxb de- 
rived with previous Suazku observations were 4-6 xlO -8 
erg cm~ 2 s _1 sr _1 (e.g. Moretti et al. 2009; Kawa- 
harada et al. 2010; Hoshino et al. 2010; Sato et al. 2012). 
The Case-GAL and Case-GAL2 gave higher CXB sur- 
face brightness by 10%. If point sources brighter than 
l.OxlO -13 erg cm" 2 s" 1 are subtracted, Scxb were mea- 
sured to be 6-7xl0 -8 erg cm -2 s" 1 sr" 1 (e.g. Bautz 
et al. 2009; Simionescu et al. 2011; Akamatsu et al. 
2011; Walker et al. 2012a,b). Using this threshold, the 
CXB surface brightness around Abcll 1835 for the Case- 
GAL+ICM is 6.8 xlO -8 erg cm -2 s _1 sr -1 which agrees 
well with the previous studies. 

3.3. X-ray Surface Brightness Profile 



Suzaku Observations of Abcll 1835 Outskirts 



5 



Case-GAL 12' -20' East 




Case-GAL 12' -20' West 



0.1 



a o.oi 



io-' 

1(H 
2 

-2 

-4 



1 2 

Energy (keV) 
Case-GAL+ICM 12" -20' East 



rt ~ b±3 -j~ r "J~-l-' J l= H 






,+ 



















0.1 



8 0.01 



10-' 

10- 1 

2 

-2 

-4 



1 2 

Energy (keV) 
Case-GAL+ICM 12" -20' West 



Energy (keV) 



Energy (keV) 



Case-GAL2 12' -20' East 



8 0.01 



10-' 

10- 4 

2 

. 
-2 
-4 



rr J=t ^^ T± ^-r-ET t >i 


^^+4-4+ 1 











0.1 



8 0.01 



10-' 

10- 4 

2 

. 
-2 
-4 



Case-GAL2 12' -20' West 



, 1 \ 


















- + 1 













Energy (keV) 



Energy (keV) 



Fig. 2. — NXB-subtracted spectra of XIS3 (black crosses) and XIS1 (red crosses) for the outermost annulus (12'— 20'). Top, middle, and 
bottom panels correspond to the Case-GAL, Case-GAL+ICM, and Case-GAL2, respectively. The ICM, CXB, LHB, MWH1, and MWH2 
emissions for the XIS3 spectra are shown in magenta, blue, orange, green, and cyan lines, respectively. Sum of the CXB, LHB, and MWH1 
emissions for the XIS3 spectra are indicated by green-cyan line. The total model spectra of XIS3 and XIS1 are shown in black and red 
lines, respectively. Left and right panels correspond to the directions in East and West, respectively. The lower panels show the residuals 
in units of a. 
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TABLE 2 

Best-fit parameters of X-ray background components 



Case 


CXB 

ScxB a 


LHB (0.1 koV) 
Norm hiiB h 


MWH1 (0.3 koV) 
Norm MWH1 h 


MWH2 (0.6 kcV) 
Afor»n MH w2 b 


12'-20' region 
X 2 /bin 


Reduced-x 2 
(xVd-o.f) 


GAL 

GAL+ICM 
GAL2 


e 7 o+0.17 
c 17 +0.26 
- CK +0.18 

5.66_ lg 


7 oq+0.63 
' -0.63 
7 54+0.66 

o , „+0.72 

S-i3_ .7i 


96+ 06 
u - MD -0.07 

77+ - 12 
u -' '-0.12 


u - iD -0.06 


694/420 
655/420 
667/420 


1.363 (2927/2148) 
1.344 (2884/2146) 
1.354 (2907/2147) 



a Estimated surface brightness of the CXB after the point source excision in unit of 10 8 erg cm 2 s 1 sr 1 (2.0-10.0 keV). 



Normalization of the apec component scaled with a factor 1/4007T assumed in the uniform-sky ARF calculation (circle radius 
r = 20'). 

Norm — 4Q Q 7r J n e nadV / [4m (1 + z) 2 D A ] x 10 -20 cm -5 arcmin -2 , where Da is the angular distance to the source. 



TABLE 3 

Estimation of CXB fluctuation for all directions 



Region 




Coverage a ' c 


SOURCEJ 1 


o-//cxB e 




(arcmin 2 ) 


(%) 


RATIO-REG (%) 


(%) 


0'-2' 


11.1 


88.1 


59.37 


23.7 


2'-4' 


33.4 


88.6 


10.90 


13.7 


4'-6' 


53.3 


84.8 


5.19 


10.8 


6'-9' 


119.2 


84.3 


4.57 


7.2 


9'-12' 


167.5 


84.6 


2.94 


6.1 


12'-20' 


590.3 


73.4 


4.01 


3.3 



a The average value of the three detectors. 

Solid angle of each observed region. 
c Fraction of each area to entire annulus. 



Fraction of the simulated cluster photons which fall in the region compared with 
the total photons generated in the entire simulated cluster. 

SOURCE .RATIO MEG = Coverage x f£ out S{r)rdr/ / °° S(r)rdr, where 5(r) rep- 
resents the assumed radial profile of Abcll 1835. We confined S(r) to a 49 49 region 
on the sky. 

c CXB fluctuation due to unresolved point sources. S c — 2xl0 -14 erg cm -2 s _1 is 
assumed for all regions. 
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In order to see how far the ICM emission of Abcll 1835 
is detected, wc derived a surface brightness profile in 
the energy band of 1-2 keV from the XIS mosaic image 
(XISO, XIS1, and XIS3 images) excluding the circular 
regions around 32 point sources (Figure 1). The left 
panel of Figure 3 shows the raw surface brightness 
profile in the 1-2 keV band (black crosses) where the 
background is included and the vignetting effect is not 
corrected. The 1-2 keV NXB profile, derived from an 
NXB mosaic image with xisnxbgen, is shown in red 
in the left panel of Figure 3. We then obtained an 
XRB mosaic image. Based on the CXB and Galactic 
(LHB+MWH1) models of Abell 1835 obtained from 
fitting for the three background cases, we simulated 
CXB and Galactic images of the four offset observations 
using xissim with exposures 10 times longer than those 
of actual observations. The 1-2 keV CXB and Galactic 
(LHB+MWH1) profiles for the Case-GAL+ICM arc 
shown in green and blue in the left panel of Figure 3, 
respectively. Since there is a contribution of luminous 
point sources outside the excluded regions, we simu- 
lated this residual-point-source signals. Details of the 
simulation are described in Appendix A. The 1-2 keV 
residual-point-source profile is shown in orange in the 
left panel of Figure 3. In the end, we obtained back- 
ground (CXB+Galactic+NXB+rcsidual-point-source) 
and background-subtracted (raw— background) profiles 
for the Case-GAL+ICM, as shown in cyan and magenta 
in the left panel of Figure 3. Here, ±10% systematic 
error for the CXB intensity is included. 

Since these images were not corrected for the vi- 
gnetting, we calculated the ratio of the background- 
subtracted surface brightness to the CXB brightness in 
the right panel of Figure 3 for the three background 
cases. As a result, this ratio decreases with radius out 
to the virial radius (r V j r — 12'0) and becomes flatter 
beyond. Beyond r v ; r , the background-subtracted signal 
in the brightness for the Case-GAL+ICM accounts for 
29% and 26% of the CXB and CXB+LHB+MWH1, re- 
spectively. It cannot be explained by the fluctuations 
of the CXB intensity (-3.3% in the 12.'0-20.'0 region 
in Table 3). An additional emission component is re- 
quired to explain the observed flux beyond r v i r . How- 
ever, it would be uncertain whether this emission come 
from the ICM (Case-GAL+ICM) or from the relatively 
hot (—0.9 keV) Galactic emission. Since Abell 1835 is 
located close to the North Polar Spur, it is clear that 
the background-subtracted signal suffers from the Galac- 
tic emission uncertainties. The ratios of the remaining 
signal to the background for the Case-GAL and Case- 
GAL2 are a factor of two smaller than those for the 
Case-GAL+ICM. From x 2 , the Case-GAL+ICM better 
reproduces the spectra beyond r v ; r . We will quantify 
the systematic errors of the ICM temperature associated 
with the background subtraction in Section 3.6 and Table 
8. 

3.4. Result of the ICM components in the All Directions 

Figure 4 shows results of the spectral fit in the particu- 
larly important regions (4.'0-6.'0, 6f0-9.'0, and 9.'0-12f0) 
toward the East and West directions (opposite azimuthal 
directions) for the Case-GAL+ICM. The best-fit param- 
eters of ICM components and x 2 values in each region, 
when parameters are tied between the four offset point- 



ings, are listed in Table 4. 

3.4.1. Temperature Profile 

Figure 5 shows radial profiles of projected temperature 
in the all directions, observed with XMM-Newton (Zhang 
et al. 2007; Snowden et al. 2008), Chandra (Bonamente 
et al. 2012), and Suzaku (this work). The temperature is 
—8 keV within 2.'0 (—480 kpc), and it gradually decreases 
toward the outskirts down to —2 keV around the virial 
radius (r v ; r — 12.'0). At a given radius, the temperatures 
for the three background cases with Suzaku agree well 
with each other. For the Case-GAL+ICM, the tempera- 
ture in the 12f0-20.'0 region joins smoothly from those of 
inner regions. We shall quantify the systematic error of 
the Suzaku ICM temperature in Section 3.6. The temper- 
ature profile derived from Suzaku agrees well with those 
from XMM-Newton within 4f0 (-0.95 Mpc) and from 
Chandra (Bonamente et al. 2012) within 7'.h (—1.8 Mpc). 
Bonamente et al. (2012) measured temperature profiles 
for Abell 1835 out to lO'O with Chandra. As their out- 
ermost temperature is in the 7.'5-10f0, we shall compare 
with our result in the West direction at the outskirts in 
subsection 3.5. 

3.4.2. Electron Density Profile 

The electron number density profile was calculated 
from the normalization parameter of apec model, defined 

as 

10~ 14 f 

NOrm = An(D A (l + z)r J ncn " dV > (2) 

where Da is the angular size distance to the source in 
units of cm, n is the electron density in units of cm~ 3 , 
and Tin is the hydrogen density in units of cnr 3 . Wc 
note that the resultant normalization using an ARF gen- 
erated by xissimarf gen needs the correction by a factor 
of SOURCE .RATIO. REG /O e (see details in Sec 5.3 of 
Ishisaki et al. 2007). The left panel of Figure 6 shows 
the radial profiles of normalization of apec model for the 
ICM component scaled with this factor. 

Each annular region, projected in the sky, includes 
emission from different densities due to integration along 
the line of sight. Assuming spherical symmetry, we 
de-convolved the normalization and calculated n c for 
each annular region from the outermost region, with the 
method described in Kriss et al. (1983). The resultant 
radial profiles of n c in the all directions with XMM- 
Newton (Zhang et al. 2007) and Suzaku (this work) are 
shown in the right panel of Figure 6 and listed in Table 
5. The deprojected electron density profile derived from 
Suzaku agrees well with that from XMM-Newton within 
5.'0 (-1.2 Mpc). The n c out to 9.'0 for the three back- 
ground cases are consistent within statistical errors with 
each other. However, n e for the Case-GAL and Case- 
GAL2 in the 9.'0-12f0 region are 26% and 17% higher, 
respectively, than that for the Case-GAL+ICM, because 
of the projection of the ICM emission outside r v j r . 

The electron density profile of Abcll 1835 for the Case- 
GAL+ICM from 6.'0 (-0.66r 200 ) to 12f0 (-1.3r 200 ) is 
well represented with a power-law model with (3 = 0.88. 
This value agrees very well with (3 = 0.89 derived for 
electron density profiles in 0.65-1.2r2oo of 31 clusters ob- 
served with ROSAT (Eckert et al. 2012). For the Case- 
GAL and Case-GAL2, the derived [3 values in the same 
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Fig. 3. — (left): Surface brightness profiles in the 1—2 keV band, using all the four observations. Point sources specified in Figure 1 
(left) arc removed, but vignetting is not corrected. Raw (background inclusive) profile (black) is shown with CXB (green), Galactic (blue), 
NXB (red), residual-point-source (orange), and background (CXB+Galactic+NXB+residual-point-source; cyan) profiles. A resultant 
background-subtracted profile (black— cyan) is shown in magenta. The error bars in these profiles are 1<t. In the error bars in the magenta 
profile, ±10% error for the CXB intensity is added in quadrature to the corresponding statistical la errors, (right): The ratio of the 
background-subtracted surface brightness (magenta of left panel) to the CXB surface brightness (green of left panel) in the 1-2 keV band 
for the three background cases. 



TABLE 4 

Best-fit parameters of ICM components in the all directions 



Region 


fcT AU a 
(keV) 


Case-GAL 
A?ormAll a ' b 


X 2 /bin c 


kT AU * 
(keV) 


Case-GAL+ICM 
A r orm,Aii a ' b 


X 2 /bin c 


fcT A11 a 
(keV) 


Case-GAL2 
AformAii a,b 


X 2 /bin c 


2'-4' 
4'-6' 
6'-9' 
9'-12' 
12'-20' 


7 e.H-0.65 
'■ oo -0.64 
4 Ql+0.99 

QQ~f"0.52 
Z.OO_q 32 
-, 79 +0.53 
x - ' z -0.15 


6.12l8;»8 x IQi 

7 9 +0.48 

— 0. 47 

1 7 4 +0.38 


682/420 
457/420 
481/420 
489/420 
694/420 


7 zz+0.64. 

1 ■°°-0.64 

c nn+ ' 95 

o.uu_ 79 

2 56+ ' 57 
z - do -0 .42 

2 00+ ' 54 

, 7 9+0.49 

1 - lz -0.36 


6.i7i8; 2 2 8 x 101 
9.i4l8:l? 

4 14+0-49 
4 - i4 -0.48 
9 90+O.4O 
z ' ZJ -0.40 

l.UO_ 27 


680/420 
457/420 
482/420 
486/420 
655/420 


7 cn+0.65 

r in + 1.07 
°- ±o -0.S6 

z - DU -0.46 
9 r l '7+0.68 
z - u ' -0.44 


6.ni8:i8 x 10 1 

8.54i8:« 

q tri+0.49 
• 3 - oi -0.47 
1 ci +0.40 
L - OX -0.38 


681/420 
458/420 
486/420 
490/420 
667/420 



T"a jj represent for temperatures 

ormalization of the apec compoi 

, m _ SOU RCE .RATIO .REG 
m _ n 

2 of the fit when p 



of the sum of all the azimuthal regions. Simile 
Lent scaled with a factor of S OU RC E.RAT I O. 
f nc n H dV I [4tt (1 + z) 2 E>\] X IO" 20 cm" 5 
ters are tied between the four offset pointings. 



way are 0.38 and 0.51, respectively, which are signifi- 
cantly smaller than those derived by Eckert et al. (2012). 
Here, r2oo is calculated using the average ICM tempera- 
ture, (kT), r 200 = 2.47/if o V(fcT)/l0keV Mpc. This re- 
lation was expected from numerical simulations for our 
cosmology (Henry et al. 2009). The average temperature 
of Abell 1835 integrated over the radial range of 70 kpc to 
r 500 (= 1.39 Mpc /ifo 1 or 5.'85) with XMM-Newtonis 7.67 
± 0.21 keV (Zhang et al. 2007), where r^oo is defined by 
weak lensing analysis (Okabe et al. 2010a). Thus, from 

the average temperature, r 2 oo = 2.16 Mpc h? or 9.'08, 
which is close to 9'29 for r2oo defined by the weak lensing 
analysis. 



3.4.3. Entropy Profile 

Figure 7 shows the entropy profiles with XMM-Newton 
and Suzaku calculated as 



rly, Norm An repn 
JJEG/O, 



TABLE 5 

Deprojected electron number densities in the all directions 



Region 


Case-GAL 


Case-GAL+ICM 


Case-GAL2 




«o,All a 


™ c ,Aii a 


™ c ,All a 



2' 


-4' 


6.64lC:i| 


X 10" 


-4 


6 - 64 -0.13 


x 10 


-4 


6 - 64 -0.13 


x 10 


-4 


1' 


6' 


i-86l8:i 4 4 


x 10" 


-4 


, ocr + 0.14 
1,8b -0.X6 


x 10 


-4 


1-86-0.14 


x 10 


-4 


6' 


-9' 


9.57+^i 


x 10" 


-5 


9.80+1? 2 


x 10 


-5 


9.56l r °? 


x 10 


-5 


9' 


-12' 


7 ni+0.83 
'' M1 -0.89 


x 10" 


-5 


6.30+i : ° 7 


x 10 


-5 


7 no+0.91 

' •• : '' 5 -o.99 


x 10 


-5 


12' 


-20' 








^• zy -0.60 


x 10 


-5 









K 



kT 

2/3 1 



(3) 



where T and n D are the temperature and deprojected 
electron density obtained above, respectively. At a given 
radius, the entropies for the three background cases 
are consistent within statistical errors with each other. 
Within ?*5oo, the derived entropy profiles with XMM- 
Newton and Suzaku follows a power-law model with a 
fixed index of 1.1, which was predicted from the accre- 
tion shock heating model (Tozzi & Norman 2001; Pon- 
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Fig. 4.— In the Case-GAL+ICM, the NXB-subtracted spectra of XIS3 (black crosses) and XIS1 (red crosses). The ICM, CXB, LHB, 
and MWH1 emissions for the XIS3 spectra are shown in magenta, blue, orange, and green lines, respectively. Sum of the CXB, LHB, and 
MWH1 emissions for the XIS3 spectra are indicated by green-cyan line. The total model spectra of XIS3 and XIS1 are shown in black 
and red lines, respectively. Left and right panels correspond to the directions in East and West, respectively. The lower panels show the 
residuals in units of a. 
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Fig. 5. — Projected temperature profiles obtained by spectral 
analyses, when all spectra of azimuthal regions at the same dis- 
tance are summed (All). Black dotted, red solid, and green 
dashed diamonds show our Suzaku results for the Case-GAL, Case- 
GAL+ICM, and Case-GAL2, respectively. The uncertainty range 
due to the ±10% variation of the CXB levels from the CXB nom- 
inal levels (red) for the Case-GAL+ICM is shown by two blue 
dashed lines. XMM-Newton results by Zhang et al. (2007) are the 
orange crosses. Weighted average of two XMM-Newton results by 
Snowden et al. (2008) are the magenta crosses. Chandra results 
by Bonamente et al. (2012) are the blue crosses. The error bars in 
XMM-Newton and Chandra results are 90% confidence level. Ver- 
tical dotted lines show rsno (5f85) and r v j r (12?0) determined by 
weak lensing analysis (Okabe et al. 2010a). 

man et al. 2003; Voit et al. 2005). In contrast, beyond 
r 500 (~6.'0), the entropy profiles for the three cases be- 
come natter in disagreement with the r 1A relationship. 
For the Case-GAL+ICM, the entropy profile is flat out 
to 20f0 (~1.7r vir ). 

3.5. Result of the Case-GAL+ICM in Each Direction 

We here investigate azimuthal variation of X-ray ob- 
servables. We fitted all the spectra simultaneously in the 
same way as for the Case-GAL+ICM (Section 3.1), but 
the temperatures and normalizations for the ICM com- 
ponent within 12.' for the four directions were indepen- 
dently determined. The derived x 2 , 2669 for 2119 degrees 
of freedom, became smaller than 2884 for 2146 degrees 
of freedom of the previous fit for all direction. Table 6 
shows the derived temperatures, normalizations, and x 2 
values for each region. Assuming spherical symmetry we 
calculated the deprojected electron number density and 
entropy profiles. The resultant radial profiles of n e in 
each direction are listed in Table 7. Figure 8 shows the 
radial profiles of the ICM temperature, scaled normal- 
ization, electron density, and entropy in each direction. 

There is an azimuthal variation of projected temper- 



atures in the outskirts of rsoo 



The best-fit 



temperature in the East region is highest, about twice of 
those in the West and North regions. The temperature of 
South region is intermediate between them. The differ- 
ences between individual values are larger than measure- 
ment uncertainties, but is not significant. Bonamente et 
al. (2012) reported with Chandra a temperature of kT 
= 1.26 ± 0.26 keV in the 7.'5-10'0 region. We note that 
they observed the western region of Abell 1835. Their 
temperature is marginally consistent with our measure- 



ments at 6f0-9.'0 (kT = 2.09l^ keV) and at 9f0-12f0 
(kT = 1.51±g;S? keV) in West direction (Figure 8 and 
Table 6). The electron number density and entropy pro- 
files in the South, West, and North directions are con- 
sistent with one another within the statistical errors in 
any annulus region. The electron number density and 
entropy in the outskirts (^500 r J" v ir) of East direc- 
tion are lower and higher than those in other directions, 
respectively. 

3.6. Systematic Errors 

We examined the effect of systematic errors on the de- 
rived spectral parameters. The level of the CXB fluctu- 
ation was scaled from the Ginga result (Hayashida et al. 
1989) as shown in Table 3. The CXB fluctuation in the 
cluster outskirts (6'0 r), where correct estimations of 
the CXB intensity are of utmost importance, is less than 
10% in the all directions from Table 3. In the 9f0-12.'0 
region, for example, the CXB fluctuation value in the all 
directions is about 6.1%. Although in the cluster center 
(r ;$ 6.'0), the CXB fluctuation is higher, the systematic 
uncertainties caused by the fluctuation became smaller 
due to much brighter ICM emission. Therefore, we as- 
sume that the upper and lower limits of the CXB system- 
atic changes in the all directions are ±10%, even consid- 
ering the uncertainty other than the systematic error due 
to the spatial variation. For the Case-GAL+ICM, we re- 
peated the spectral fit for all directions in the same way 
but fixed the CXB intensity at the upper and lower 10% 
from their nominal levels. Table 8 shows the changes of 
ICM properties in the 9'0-12f0 region and the reduced- 
X 2 in each variation. In the 9.'0-12f0 region, the effects of 
±10% error for the CXB intensity, on the temperature, 
electron number density, and entropy are 20%-30%, less 
than 5%, and 20%-30%, respectively. In each direction, 
the CXB fluctuation for each annular region should be 
a factor of two higher than that of for all directions. In 
the 9'0— 12'0 region, for example, the CXB fluctuation 
values in the East, South, West, and North directions 
are about 11.5%, 12.6%, 11.2%, and 10.6%, respectively. 
Therefore, we assume that the upper and lower limits of 
the CXB systematic changes in each direction are ±15%. 
In each direction, the change of entropy in the 9f0-12.'0 
region is 20%-40% by ±15% error for the CXB intensity, 
except for the CXB+15% in the West direction. This 
systematic error is comparable to the statistical error. 
For the CXB+15% in the West direction, the tempera- 
ture and entropy are lower by a factor of ^2 and ^3, 
respectively. 

The changes of the intensities of the Galactic emissions 
(LHB, MWH1), by systematic changes in the CXB lev- 
els to ±10% and ±15% from their nominal levels, arc 
less than 10%. Then, we repeated the spectral fit by fix- 
ing the Galactic (LHB, MWH1) intensity at the upper 
and lower 10% from their nominal levels, for the Case- 
GAL+ICM, and by fixing the metal abundance of the 
ICM at 0.1 and 0.3, instead of 0.2. We also show these 
results in Table 8. The systematic errors due to these 
effects are less than that for the CXB in the 9f0-12.'0 
region. 

Since the spectra with azimuthal variations in temper- 
ature are focused to fit with a single temperature model, 
we obtain relatively large reduced-% 2 , over 1.3 (see Table 
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Fig. 6. — (left): The radial profiles of normalization of apec model obtained by spectral analyses scaled with a factor of 
SOU RC EJiATIO-REG /Q e , when all spectra of azimuthal regions at the same distance are summed (All). Dotted (black), solid (red), 
and dashed (green) diamonds show our Suzaku results for the Case-GAL, Case-GAL+ICM, and Case-GAL2, respectively. The uncertainty 
range due to the ±10% variation of the CXB levels from the CXB nominal levels for the Case-GAL+ICM is shown by two dashed (blue) 
lines. Vertical dotted lines show rsoo (5'85) and r v ; r (12'0). (right): The same as left panel, but for deprojected electron number density 
profiles. See the text for the detailed method of derivation (Section 3.4.2). XMM-Newton results by Zhang ct al. (2007) are the crosses 
(orange) . 



TABLE 6 

Best-fit parameters of ICM components in each direction {Case-GAL+ICM) 



Region 


fcT East a 


A r orm East a ' b 


fcT Sout h a Norm South a ' b 


fcT Wcat a 


7Vor-mwo S t a ' b fcT Nort h a 7Vorm Nortn a ' b X 2 /bin c 




(keV) 




(keV) 


(keV) 


(keV) 



0' 


-2' 










o ocr+0.63 
°- zo -0.53 


1.35i°-gt x 10 3 






125/105 


2' 


1' 


7.74+J-|° 


4.83±° g x 10 1 


7.24l?i° 


4.551°;™ x 10 1 


D.DO_q 7g 


5.29l° :3 ; x 10 1 


8.6o±; : g 


8.18±°; 3 33 x 10 1 


495/420 


4' 


6' 


6.40t?il 


9.13±M? 


o q 5 +1.74 




d - du -1.32 


»-os±i:|g 


3.98+i-n 




452/420 


6' 


-9' 


3 - 43 ±i'.17 


o o +0.89 


z -'°-0.64 


c 11+0.98 
0,1 -0.94 


o nq+ - 68 
z -uy_ 4 7 


o 74 +0.85 
°- '^-0.80 


2 zn+ a7s 
z - 4J -0.58 




471/420 


9' 


-12' 


4.01_2.34 d 


-, 1O+0.77 
i - L8 -0.6S 


9 S1+ 2 ' 11 
z - 81 -0.91 


, no +0.92 
o.UZ_ Q g5 


i ci +0.49 
1 - O1 -0.27 


9 .. . +0.69 
z - o± -0.67 


l-84_ 32 




471/420 



fcT East , fcTg outh , fcT'westi and fcT Nortn represent for temperatures of East, South, West, and North, respectively. Similarly, Norm-^ 
N OT ' m ]SJoi : th rc P reselr t f° r normalization. 

b Normalization of the apec component scaled with a factor of SOU RC E.RAT I O.REG / 'Ct c 

Norm = SOURCE.RATIO.REG J^^rfy/ [4^ (1 + ,) 2 D 2 ] X 10" 20 cm" 5 arcmin" 2 , where Z? A is the angular distance to the source 
C of the fit when parameters are separately determined for the four offset pointings. 

Wo could not constrain the 90% confidence level upper limit of the temperature because of the influence of bright point sources (sec Appc 



TABLE 7 

Deprojected electron number densities in each direction (Case-GAL+ICM) 



Region n CiEas t a ™ ,South a n c ,Wost a n e , North 8 



0' 


-2' 














4 9Q+ - 07 
*- zy -0.07 


x 10" 


3 








2' 


4' 


r- -ro+0.43 
°- <°-0A8 


x 10" 


-4 


r, ^4+0.52 
°-°^-0.49 


x 10" 


-4 


°- u ' -0.21 


x 10" 


-4 


7 7 o+0.18 
'•'O-0.18 


x 10" 


-4 


■1' 


6' 


9 i 9 +0.21 
z - iz -0.22 


x 10" 


-4 


, 7t -+0.30 


x 10" 


-4 


1 qi+0.22 
1 - al -0.23 


x 10" 


4 


i S9 +0.31 
± -° z -0.36 


x 10" 


-4 


6' 


-9' 








l.UO_ 22 


x 10" 


-4 


8.7o±J:i; 


x 10" 


-5 


1 02+ - 18 
i - uz -0.21 


x 10" 


-4 


9' 


-12' 










x 10" 


-5 


6.541J;™ 


x 10" 


5 


6-93ii;U 


x 10" 


-5 


6' 


-12' 


^ c;«+ 1 - 15 
°-°°-1.39 


x 10" 


-5 





















a ™c,East> ™c, South' n c,Wcst> and n e, North represent for electron number densities of East, South, West, 
and North, respectively, in units of cm -3 
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Fig. 7. — The radial profiles of entropy in the all directions ob- 
tained by calculating K = kT/n^ 3 from profiles in Figure 5 and 
right panel of Figure 6. Dotted (black), solid (red), and dashed 
(green) diamonds show our Suzaku results for the Case-GAL, Case- 
GAL+ICM, and Case-GAL2, respectively. The uncertainty range 
due to the ±10% variation of the CXB levels from the CXB nom- 
inal levels for the Case-GAL+ICM is shown by two dashed (blue) 
lines. Vertical dotted lines show rsoo (5.'85) and r vir (12.'0). XMM- 
Newton results by Zhang et al. (2007) are the crosses (orange). 
The solid (orange) line extrapolates the XMM-Newton data with 
a power-law formula fit for the data beyond 70 kpc, which has a 
fixed index of 1.1. 

2 and Table 8). A part of the large-x 2 is due to the az- 
imuthal variation in the temperature and normalization 
at the 2f0-4.'0 region: when we set these parameters in 
the four azimuthal directions to be free, the x 2 for this 
annular region improved to 495 (Table 6) from 680 (Ta- 
ble 4) for the all directions. Furthermore, we estimate 
changes of x 2 by adding 8% systematic error in the spec- 
tral data points for uncertainty analysis. The minimum 
reduced-x 2 (x 2 /d.o.f) improved to 1.081 (2321/2146) in 
the all directions for the Case-GAL+ICM. There was no 
significant effect on the results of the 6'0-9f0 and 9.'0- 
12'.0 regions. In particular, the minimum \ 2 f°r the 9f 0— 
12'0 region significantly improved from 486 to 412 for 
420 bins. Accordingly, the ICM temperature and nor- 
malization of this region slightly changed from 2.00±g;|| 
kcV to 2.03+S-Sn keV and from 2.23±0.40 x 10~ 20 cm" 5 



-0.40 

'- 2 to 2.19+S-fl x lO" 20 



cm 5 arcmin 2 , respec- 



aremm 
tively. 

We next discuss the influence of Suzaku's PSF. We ex- 
amined how many photons accumulated in the six an- 
nular regions actually came from somewhere else on the 
sky because of the extended telescope PSF following a 
procedure described by Sato et al. (2007). Table 9 shows 
the contribution from each sky region for 0.5-2 keV en- 
ergy range. Here, we averaged the values of the three 
detectors. Although the effect of leakage from the cen- 
ter < 2f0 to the 2?0-4.'0 and 4f0-6f0 regions is severe, 
the derived temperature and normalization in this re- 
gion agree well with those derived from XMM-Newton 
(Zhang et al. 2007; Snowden et al. 2008) and Chandra 
(Bonamentc et al. 2012). In the cluster outskirts, out- 
side 6'0, the contributions of the leakage from the central 
4f0 region are less than 20%. We fitted the spectra for 



the 9.'0-12f0 region, including the stray light component 
from the bright central region. Then, we got 8% lower 
ICM temperature and 9% smaller normalization. Even if 
the actual effect of the stray light was twice the current 
calibrations, the changes would be —14% and —16% for 
the temperature and normalization, respectively. These 
differences are significantly smaller than the present sta- 
tistical errors. Therefore, the temperature changes by 
the PSF correction should be mostly small as in previ- 
ous Suzaku observations of cluster outskirts (George et 
al. 2009; Reiprich et al. 2009). 

Although the temperature profile is a projected one 
obtained by the two-dimensional spectral analysis, the 
results of the deprojection fitting do not show any signif- 
icant difference from the non-deprojection fitting (Bautz 
et al. 2009; Akamatsu et al. 2011; Walker et al. 2012a). 

4. DISCUSSION 

Suzaku performed four pointings, deep observations of 
Abell 1835 toward four azimuthal directions (Figure 1). 
Faint X-ray emission from the ICM out to the virial ra- 
dius was detected, enabled us to measure radial pro- 
files of gas temperature, electron density, and entropy 
The X-ray observables outside the virial radius, for the 
Case-GAL+ICM, are in good agreement with extrapo- 
lated from the inside, although we cannot distinguish 
whether the emission comes from the ICM or from rel- 
atively hot Galactic emission. We here discuss cluster 
thermal properties within the virial radius, incorporating 
weak-lensing mass from Subaru/Suprime-Cam, Sunyaev- 
Zel'dovich (SZ) effect flux from Planck satellite, and pho- 
tometric data from SDSS. 

4.1. A joint X-ray and Lensing Analysis 
4.1.1. Mass Estimation 

We here estimate the hydrostatic equilibrium (H.E.) 
masses for three background models using azimuthally 
averaged temperature and electron density profiles and 
compare weak-lensing mass (Okabe et al. 2010a). We 
first fitted the temperature profiles from XMM-Newton 
(from 0'.5 to 5.'0; Zhang et al. 2007) and Suzaku (from 2.'0 
to 20.'0) results with the (3 model with constant, and the 
electron density profiles from XMM-Newton (from O.'O 
to 3.'0) and Suzaku (from 2.'0 to 20f0) results with the 
double-/? + power-law model. The results for the Case- 
GAL+ICM are plotted in Figure 9. The hydrostatic 
mass, Mh.e.(< r), within the three-dimensional radius 
r, assuming spherically symmetric hydrostatic equilib- 
rium, is calculated from the parametric temperature and 
electron density profiles with the following formula (Fab- 
ricant et al. 1980), 



M H . E .(< r) = - 



kT(r)r ( <ilnp g (r) dhiT(r] 



(irripG \ dhir 



dhir 



,(4) 



where G is the gravitational constant, k is the Boltzmann 
constant, p g is the gas density, kT is the temperature, 
and [im v is the mean particle mass of the gas (the mean 
molecular weight is fjb = 0.62). We also calculated gas 
mass M gas by integrating the electron density profile. 

Weak-lensing analysis with the Subaru/Suprimc-Cam 
is described in detail in Okabe et al. (2010a). Weak- 
lensing mass, Mi ons , is derived by fitting Navarro Frcnk 
& White (NFW) model (Navarro, Frenk, & White 1997) 
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Fig. 8. — Comparison of temperature, normalization scaled with the factor, deprojected electron number density, and entropy profiles in 
the East (black), South (red), West (green), and North (blue) observations for the Case-GAL+ICM. All directions are magenta crosses. 
XMM-Newton results by Zhang et al. (2007) are the orange crosses. 



TABLE 8 

The changes of ICM profiles in the 9.'0-12.'0 region and the x 2 in each variation by 

SYSTEMATIC CHANGES IN THE ALL DIRECTIONS (Case-GAL+ICM) 



Variation Temperature Density Entropy Reduced-x 2 (x 2 /d.o.f) Ax 2 " 



CXB+10% 


-14% 


+2% 


CXB-10% 


+29% 


+1% 


Galactic+10% 


+8% 


+1% 


Galactic- 10% 


-13% 


-0% 


abundance 0.1 (fixed) 


-14% 


+6% 


abundance 0.3 (fixed) 


+5% 


-3% 


The difference in % 2 from that fc 


r the Ca. 


e-GAL+ICM 



to the tangential lensing distortion profile in the range 
of lf0-18.'0. The NFW profile is well defined form for 
the spherically averaged density profile of dark matter 
halos with high-resolution numerical simulations. They 
are parametrized by the virial mass (M v ; r ) and the halo 
concentration parameter (c v ; r ). 

Figure 10 shows radial profiles for the H.E. masses 
(Mh.e.) and gas one (M gas ) using the three background 
cases (Case-GAL+ICM, Case-GAL2, and Case-GAL), 
and weak-lensing mass. The mass (Mh.e. and M gas ) pro- 



-15% 


1.348 (2895/2147) 


+11 


+29% 


1.347 (2891/2147) 


+7 


+8% 


1.348 (2896/2148) 


+12 


-13% 


1.347 (2894/2148) 


+10 


-17% 


1.347 (2890/2146) 


+6 


+7% 


1.342 (2880/2146) 


-4 



files out to 9'.0 for the three cases are consistent within 
statistical errors with each other. Mh.e. outside 8.'0 
unphysically decreases with increasing the radius. This 
means that the hydrostatic equilibrium we assumed is 
inadequate to describe the ICM in the outskirts. The 
Mh.e. agrees with Mi ens within l.'0-5.'0 from the cluster 
center, but there is a significant difference in the regions 
of r < l.'O and r > r^oo ~ 5'85. The H.E. mass is lower 
than weak-lensing one outside rsoo because of the break- 
down of hydrostatic equilibrium assumption. The weak- 
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TABLE 9 

Estimated fractions of the ICM photons accumulated in 

DETECTOR REGIONS COMING FROM EACH SKY REGION FOR ALL 
DIRECTIONS (0.5-2 KEV) 



Detector region 


0'-2' 


2'-4' 


Sky region 
4'-6' 6'-9' 


9'-12' 


12'-20' 


0'-2' 


96.7% 


3.1% 


0.1% 


0.0% 


0.0% 


0.0% 


2'-4' 


57.2% 


37.2% 


5.2% 


0.4% 


0.0% 


0.0% 


4'-6' 


20.5% 


23.4% 


46.0% 


9.6% 


0.4% 


0.1% 


6'-9' 


15.7% 


5.5% 


15.9% 


55.7% 


6.8% 


0.5% 


9'-12' 


10.5% 


2.8% 


2.3% 


17.9% 


58.0% 


8.4% 


12'-20' 


2.0% 


1.0% 


1.2% 


2.4% 


11.6% 


80.6% 



lensing mass inside l'.O, extrapolated from the best-fit 
model obtained by fitting in the region of l.'O < r < 18.'0, 
is significantly lower than H.E. mass. As weak-limits of 
lensing distortion is breakdown in this region, the strong 
lensing method plays an important role to reconstruct 
mass distribution. A further study of a strong- and weak- 
lensing joint analysis, as demonstrated by Broadhurst et 
al. (2005), would be powerful to obtain lensing mass dis- 
tribution for the entire radial range. 

4.1.2. Gas mass fraction 

We derived the cumulative gas mass fraction within 
the three-dimensional radius r. defined as 



/gas(< r) 



M gas (< r) 
« rY 



AL 



total 



(5) 



where M gas (< r) and M to tai(< r ) are the gas mass and 
the gravitational mass (hydrostatic mass Mh.e. or lens- 
ing mass Miens), respectively. 

The gas mass fraction, fg^iK r) = M sas /M\ ens , 
combined with complementary Suzaku X-ray and 
Subaru/Suprime-Cam weak- lensing data set, is shown in 
Figure 11. As lensing mass docs not require an assump- 
tion of hydrostatic equilibrium, a comparison of results 
derived solely by X-ray data with by joint analysis al- 
lows us to understand the ICM states and the system- 
atic measurement bias. We found no significant differ- 
ence between the three background cases. The gas mass 
fraction in the range of r25oo ~ t ~ r v \ x is approximately 
constant, accounting for ^90% of cosmic mean baryon 
fraction, Q,\,/£l m , derived from seven-year data of Wilkin- 
son microwave anisotropy probe (WMAP 7; Komatsu et 
al. 2011). It is in good agreement with recent numerical 
simulations (Kravtsov et al. 2005; Nagai et al. 2007). It 
is also consistent with J7b/f2 m at r v - n - within large er- 
ror. There is no significant radial dependence of gas 
mass fraction, to the contrary of recent statistical studies 
using lensing and X-ray data set (Mahdavi et al. 2008; 
Zhang et al. 2010). The best-fit value reaches f2b/f2 m 
at r ~ l.lr v i r (r = 13.'0). The gas mass fraction in the 
central region (r < l'.O) increases to the cluster center, 
conflicting with numerical simulations (Kravtsov et al. 
2005; Nagai et al. 2007) showing that the cooling pro- 
cess increases the stellar fraction and, correspondingly, 
decreases the gas mass fraction. It implies that the lens- 
ing mass is underestimated in the central region in which 
strong lensing data is essential to reconstruct mass dis- 
tribution. Future study of joint strong and weak lensing 
analysis will allow a detailed examination in the cluster 
center. 



The gas mass fraction derived by X-ray data, /g as (< 
r) = -A/gas/A/H.E., increases monotonically with increas- 
ing radius. The /g£' out to the intermediate ra- 
dius (r < 0.5r5oo) agrees with numerical simulations 

(Kravtsov et al. 2005; Nagai et al. 2007) and f^ s) - It, 
however, significantly exceeds tt-b/Q m beyond rsoo, as 
found by the Suzaku observations of the north-west di- 
rection of the Perseus cluster (Simionescu et al. 2011). 
They interpreted the significant excess as the presence of 
gas-clumping, mainly in the cluster outskirts (Nagai & 
Lau 2011). However, our gas mass fraction estimated 
by gas and lensing masses in the range of f5oo ~ r ~ r v - n -, 
does not exceed the cosmic mean baryon fraction from 
WMAP, which indicates that the gas-clumsiness effect 
mentioned by Simionescu et al. (2011) is less significant. 
The gas mass fraction, along with the low temperature 
and entropy in the cluster outskirts, support that the 
underestimate in H.E. mass is due to the breakdown of 
hydrostatic equilibrium in cluster outskirts, rather than 
the gas-clumpiness effect. To balance fully the gravity of 
lensing mass, we need additional pressure supports such 
as turbulences and bulk motions caused by infalling mat- 
ter at the cluster outskirts. The other possibility is devi- 
ations in the electron and ion temperatures at the cluster 
outskirts (e.g. Takizawa 1999; Hoshino et al. 2010; Aka- 
matsu et al. 2011). This possibility is discussed in Section 
4.5. 

4.2. Comparison of Temperature Profile with Other 
Systems 

The temperature was measured out to the virial radius 
r v ; r ^2.9 Mpc, which shows that it decreases from about 
8 keV around the center to about 2 kcV at r v ; r . Figure 
12 compares azimuthally averaged temperature for the 
Case-GAL+ICM with the fitting function for the out- 
skirts temperature profile derived from numerical simu- 
lations (Burns et al. 2010). The fitting function, 



T 



= A 



1 + 

V r 200 / 



13 



(6) 



is scaled with average temperature and r2oo, and they 
obtained the best-fit values of A = 1.74 ± 0.03, B = 
0.64 ± 0.10, and /3 = -3.2 ± 0.4. The temperatures out- 
side 0. 3r2oo agree with that of Burns et al. (2010) within 
lcr error range (dashed lines) and those of other clus- 
ters observed by Suzaku (Burns et al. 2010; Akamatsu 
et al. 2011). It implies that low temperature (^0.2(fcT)) 
in cluster outskirts is a common feature. Burns et al. 
(2010) have shown that the kinetic energy of bulk and 
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Fig. 9. — (left): The radial profile for azimuthally averaged temperature from XMM-Newton (from 0.5.'0 to 5.'0; Zhang et al. 2007) and 
Suzaku (from 2'0 to 20'.0) for the Case-GAL+ICM. The black solid line and gray regions are the best-fits and the 68% CL uncertainty 
error ranges, respectively, (right): The same as left panel, but for the electron density from XMM-Newton (from O.'O to 3.'0) and Suzaku 
(from 2'.0 to 20f0). 
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Fig. 10. — Comparison of the hydrostatic (M tota i) and weak- 
lensing (Mi cns ; Okabe et al. 2010a) masses, with the gas mass 
(M gas ). The upper black, red, and green solid lines represent the 
total mass M tota i for the Case-GAL, Case-GAL+ICM, and Case- 
GAL?,, respectively. The lower black, red, and green solid lines 
represent the same as upper lines, but for the gas mass, M gas . 
Blue solid line represents the weak lensing mass. The solid and 
dashed lines are the best-fit values and the 68% CL uncertainty 
errors, respectively. 

turbulent motions is as high as 1.5 times of the thermal 
energy at the virial radius. As reported by recent nu- 
merical simulations (e.g. Nagai et al. 2007; Piffaretti 
& Valdarnini 2008; Jeltema et al. 2008; Lau et al. 2009; 
Fang et al. 2009; Vazza et al. 2009; Burns et al. 2010), the 
kinetic pressure in the ICM would contribute to balance 
the total mass. 

4.3. Comparison of Entropy Profile with Other Systems 

Sato et al. (2012) compared the entropy profiles with 
several clusters of galaxies observed with Suzaku. When 
scaled with average ICM temperature, the entropy pro- 
files for clusters with ICM temperatures above 3 keV are 
universal irrespective of the ICM temperature. Walker 
et al. (2012c) confirmed that the function represents well 



Fig. 11. — Cumulative gas mass fraction, / ga s(< r), averaged 
over all azimuthal directions. The solid and dashed lines represent 

the results obtained with the lensing-based total mass, /|La ns '(< 
r) = Mgas/Mi ens , and with the H.E. total mass, /gS' E ''(< r) = 
A^gas/Afe.E.j respectively. Black, red, and green solid lines rep- 
resent the gas mass fraction for the Case-GAL, Case-GAL+ICM, 
and Case-GAL2, respectively. The inner and outer two lines are 
the best-fit values and the 68% CL uncertainty errors, respectively. 
Horizontal solid two lines show the error range of the cosmic mean 
baryon fraction, fl b /Q m from WMAP (Komatsu et al. 2011). The 
^2500 (= 0.54 Mpc h^ or 2.'29) is defined by weak lensing analysis 
(Okabe et al. 2010a). 

ICM entropy profiles of several clusters of cluster, in 
agreement with the entropy profiles obtained by combin- 
ing Plank, ROSAT, and XMM-Newton. The entropy pro- 
file of Abell 1835 also becomes flat beyond 0.5r2oo, con- 
trary to the r relationship expected from the accretion 
shock heating model (Tozzi & Norman 2001; Ponman et 
al. 2003; Voit et al. 2005). One possible explanation for 
the low entropy profiles at cluster outskirts is that ki- 
netic energy accounts for some fraction of energy budget 
to balance fully the gravity (Bautz et al. 2009; George 
et al. 2009; Kawaharada et al. 2010; Sato et al. 2012). 
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Fig. 12. — Comparison of the temperature profiles for the Case- 
GAL+ICM and the function by Burns et al. (2010). The profiles 
are scaled with the average ICM temperature, (kT), and T200 value 
derived from Henry et al. (2009) for fair comparison with previous 
studies. Circles (magenta) and boxes (orange) show our Suzaku 
result and XMM-Newton result (Zhang et al. 2007), respectively. 
Errors are the 68% CL uncertainty. The dotted line shows the 
simulation result of Burns et al. (2010). Two dashed lines show 
the standard deviation. 

The flattening of the entropy profile supports that the 
discrepancy between the Mh.e and Mi ens beyond 0.5r v ; r 
is caused by the deviation from the hydrostatic equilib- 
rium. 

4.4. Comparison with Planck stacked SZ pressure profile 

X-ray emission is proportional to the square of the 
gas density integrated along the line of sight, thus it 
is powerful for the denser region of the hot gas. The 
thermal Sunyaev-Zel'dovich (SZ) effect (Sunyaev & Zel- 
dovich 1972) is proportional to the thermal gas pressure 
integrated along the line of sight. The SZ observation 
therefore makes a powerful diagnostic of the less dense 
gas, like in cluster outskirts. As the sensitivity of the 
SZE to gas-clumping (Nagai <fc Lau 2011) is a func- 
tion of the pressure differential of the clumps with the 
surroundings, X-ray and SZE obscrvablcs arc thus com- 
plementary and allow us to further constrain the physics 
of the ICM. Planck is the only SZ experiment with a full 
sky coverage, able to map even nearby clusters to their 
outermost radii and offering the possibility of an in-depth 
statistical study through the combination of many obser- 
vations (Planck Collaboration et al. 2012). We compare 
pressure profile, P = kTn e , using Suzaku observed pro- 
jected temperature and deprojected electron density for 
the Case-GAL+ICM with that derived by stacked SZ 
flux for 62 nearby massive clusters from the Planck sur- 
vey (Planck Collaboration et al. 2012). The results are 
plotted in Figure 13. Following Planck Collaboration et 
al. (2012), we normalized the pressure at ?'5oo = 1-39 
Mpc /lyp 1 measured by weak lensing analysis (Okabe et 
al. 2010a). The outskirts thermal pressure measured by 
Planck agrees with our Suzaku X-ray measurement from 
2'.0 (-0.3r 500 ) to 12.'0 (~2r 50 q ~ r vir ), albeit different 
sensitivities of gas density, which indicates the reliabil- 
ity of Suzaku measurements. Walker et al. (2012c) have 
shown this agreements for other clusters observed by 
Suzaku. The consistency between independent obscrv- 



FlG. 13. — Comparison of the pressure profiles of Abcll 1835 and 
the Planck stacked SZ result of Planck Collaboration et al. (2012). 
The profiles are scaled with P500 and /(M) (see Eq.10 and Section 
6.2. in Planck Collaboration et al. 2012). and are azimuthally 
averaged. The (magenta) crosses and lines show our Suzaku result 
from the observed projected temperature and deprojected electron 
density. The error bars and the inner and outer two (magenta) lines 
are statistical lcr errors. Dashed (red), dotted (green), and chain 
(blue) lines show the best fit NFW model combined XMM-Newton 
and Planck pressure profile with all clusters, cool-core clusters, 
and non cool-core clusters, respectively (Planck Collaboration et 
al. 2012). 

ables implies that there is no strong need to invoke gas- 
clumping effect discussed by Simionescu et al. (2011). 

4.5. Comparison of ICM Profiles with Abell 1689 

Gravitational lensing masses are available for two clus- 
ters, Abell 1835 (Okabe et al. 2010a) and Abell 1689 
(Umetsu & Broadhurst 2008), among those of which X- 
ray emissions entirely within the virial radius are de- 
tected with Suzaku. A gravitational lensing on back- 
ground galaxies enables us to directly reconstruct the 
mass distribution without resting on any assumptions 
on the relation between dark matter and baryon distri- 
butions. It is important for understanding gas properties 
to compare X-ray observables based on lensing proper- 
ties. Lensing distortion profiles for two clusters are well 
expressed by the universal NFW model (Navarro, Frenk, 
& White 1997) but not well by an singular isothermal 
sphere (SIS) model. The virial masses and halo con- 
centrations for the NFW model are M vir = 1.37^:29 x 
IO^Mq/i" 1 and c vir = 3.35tg;?j) for Abell 1835 (Ok- 
abe et al. 2010a) and M vir = 1.47±S;|1 x IO^Mq/i" 1 
and c vir = 12.7 ± 2.9 for Abcll 1689 (Umetsu & Broad- 
hurst 2008; Kawaharada et al. 2010), respectively. The 
virial masses for two clusters are similar, whereas the 
concentration parameter for Abell 1835 is lower than 
Abell 1689. Indeed, Einstein radius determined by strong 
lensing analysis (Richard et al. 2010) for Abell 1835 
(3075) is smaller than that of Abell 1689 (4771) in a case 
of source redshift at z s = 2. It is well established from 
CDM numerical simulations that the halo concentration 
is correlated with the halo formation epoch (e.g. Bullock 
et al. 2001). This is because the central mass density 
for clusters, corresponding to the concentration, is corre- 
lated with those for their progenitors. Therefore, among 
clusters with similar mass, the age of clusters with higher 
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concentration is likely to be longer (Fujita & Takahara 
1999). 

If the low thermal pressure and entropy in cluster 
outskirts discovered by Suzaku (e.g. Kawaharada et al. 
2010; Sato et al. 2012; Walker et al. 2012c) are explained 
by a difference between electron and ion temperatures 
(e.g. Takizawa 1999), it is statistically expected that the 
thermalization between electrons and ions through the 
Coulomb interaction goes on for clusters with high con- 
centration. This is because the thermal equilibration 
time (< 1 Gyr) in cluster outskirts is shorter than their 
ages, although Wong & Sarazin (2009) have shown that 
electron and ion temperatures differ by less than a per- 
cent within the r v ; r . As both the virial mass and redshift 
for two clusters are similar, they are a good sample to 
compare the ICM properties and investigate a depen- 
dence of halo concentration. 

Figure 14 shows a comparison of temperature, electron 
number density, and entropy profiles in the all directions 
for Abell 1835 and Abell 1689, where the temperature 
and entropy are scaled with the average ICM temper- 
ature, (kT). The radius is also normalized by r v ; r de- 
termined by lensing analyses. They all are consistent 
within statistical errors. It implies that the temperature 
difference is not critical role to explain the observed ICM 
properties in cluster outskirts. As we don't know details 
of formation history, a statistical approach is essential to 
conclude it. 

4.6. Correlation between outskirts temperature and 
large-scale structure 

We here discuss the correlation between temperature 
anisotropy in cluster outskirts and large-scale structure 
surrounding the cluster. Kawaharada et al. (2010) have 
discovered in Abell 1689 that the anisotropic distribu- 
tions of gas temperature and entropy in the outskirts 



(»"500 



r ) are clearly correlated with the large- 



scale structure. Similarity, we found anisotropic distribu- 
tions of temperature and entropy in the outskirts, though 
their significance levels are low. 

We first make a visualization of galaxy distribution 
using photometric redshifts from the SDSS DR7 catalog 
(Abazajian et al. 2009), 



n(0) = J2 / dO' Pi (z,e')w(\e-e'\), 



(7) 



where p%(z, 6') is a photometric redshift probability dis- 
tribution for i-th galaxy and W(9) is a weight function. 
The probability function is computed from the best-fit 
z pho t and error a z , p(z, 6') = Aexp[-(x - z phot ) 2 /2/a% 
ignoring a secondary solution of photometric redshift, 
where A is the normalization by making the integration 
from z — to z — oo equal to unity. We here apply Gaus- 
sian smoothing function W(6) = l/(irQ 2 g ) cxp(-|0| 2 /0 2 ) 
with the angular smoothing scale 9 g . We used a half 
of virial radius g = lO.'O (FWHM = W.7) for recon- 
structing the two-dimensional map of galaxies. We select 
bright galaxies with the magnitude r' < 22 in a photo- 
metric redshift slice of, \z — z c \ < Sz = cr„ imax (l+2 c )/c ~ 
0.0125, where z c is the cluster redshift, z is a photometric 
redshift, cr V)lnax = 3000 km/s, and c is the light velocity. 
The resultant map is shown in Figure 15. The white 
circle and boxes represent the virial radius obtained by 
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Fig. 14. — Comparison of temperature (top), deprojected elec- 
tron number density (middle), and entropy profiles (bottom) be- 
tween Abell 1689 (Kawaharada et al. 2010) and Abell 1835 (this 
work, the Case-GAL+ICM). The temperature and entropy pro- 
files are scaled with the average ICM temperature, (kT), and az- 
imuthally averaged. The scaling radius r v ; r are determined by 
lensing analysis. Dashed (black) and solid (red) diamonds show 
the profiles of Abell 1689 and Abell 1835, respectively. XMM- 
Newton results of Abell 1835 by Zhang et al. (2007) are the crosses 
(orange) . 

weak-lensing analysis and the FoVs of the XIS point- 
ings, respectively. A filamentary overdensity region out- 
side the virial radius is apparently found in East and 
South regions of the XIS pointings. The broad filamen- 
tary structure is elongated out to ^4r V j r in the direction 
of South region. The galaxy number around the virial 
radius in East region is higher than that in South re- 
gion. The best-fit temperatures in an outermost region 
for these two regions are higher than those for the other 
region in Table 6. The North and West regions with low 
temperature contact with low density void environments. 
The correlation between the temperature in cluster out- 
skirts and the large-scale structure outside Abell 1835 is 
consistent with the result of Abell 1689 (Kawaharada et 
al. 2010). 

In order to investigate more precisely this correlation, 
we compare the outskirts temperatures for these two 
clusters with the number density contrast, 5 = nj (n) — 1, 
of large-scale structure, without any smoothing proce- 
dures applied in the map-making. The number densi- 
ties n in the azimuthal angles of the FoVs of Suzaku are 
computed by eq. 7 with W{&) = 1 and the area nor- 
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malization. Here, (n) is azimuthally averaged number 
density. We measure the density contrast in the annu- 
lus regions (l-2)r v i r and (2-4)r v i r centering the brightest 
cluster galaxy (BCG). The photometric redshift slices 
are calculated by a v max = 3000 km/s. The standard 
errors are estimated by the bootstrap method. We mea- 
sure the deviation of temperatures, T/T\\\ — 1, with the 
temperature T a n in the all directions. For a consistency 
with X-ray analysis of Abell 1689 (Kawaharada et al. 
2010), we use for Abell 1835 the temperatures measured 
in 6f < r < 12'.Q which is close to rsoo ~ r r v - u . There 
is an apparent correlation between two quantities (Fig- 
ure 16). In order to quantify this, we computed Spear- 
man's rank correlation coefficient, r s . The errors are es- 
timated by 10 4 Monte Carlo redistributions, taking into 
account the uncertainties. Full dataset for two clusters 
gives r s = 0.762 ± 0.133 for (l-2)r vir (Table 10), and the 
probability of a null hypothesis, P = 0.028t°^tl, that 
there is no relationship between two data sets. It indi- 
cates that the correlation is not occurred by chance with 
more than 90% confident. The result does not change by 
choosing luminosity- weighted center or number- weighted 
center within the virial radius. We also confirmed the 
same results using the temperature and entropy in the 
outermost regions (Table 6). The Spearman's coefficient 
using the regions of (2-4)r v i r is smaller and the proba- 
bility of a null hypothesis is more than significance level 
of 10%. The correlation with galaxy distribution in the 
large-scale structure at distance of 10 Mpc is not statisti- 
cally strong. We also investigated Spearman's coefficient 
for each cluster and could not rule out a null hypothesis, 
as expected from a low significance level of temperature 
anisotropic distribution. 

Combining the results of two clusters, we found that 
the anisotropic temperature distribution in the cluster 
outskirts is significantly associated with contacting re- 
gions of large-scale structure environment. Further study 
with larger sample is of vital importance to obtain more 
robust result. 

5. SUMMARY 

We observed Abell 1835 (temperature ~ 8 keV) with 
Suzaku and detected the ICM emission out to the virial 
radius, r v i r (~2.9 Mpc or 12.'0). Surface brightness pro- 
files and results of the spectral fit need an emission com- 
ponent in addition to those of the CXB and Galactic 
(LXB+MWH1) beyond r v i r , although we are not able to 
distinguish whether this emission comes from the ICM or 
from relatively hot Galactic emission. We derived radial 
profiles of temperature, electron density, and entropy, at- 
tempted to evaluate the cluster properties in the cluster 
outskirts, and discussed their implications. Our conclu- 
sions in this work are summarized as below. 

• The temperature gradually decreases with radius 
from ~8 keV in the inner region to ~2 keV at r v ; r . 

• The electron density profile continuously steepens 
with radius out to r vlI . The slope of the electron 
density of the Abell 1835 agrees well with those of 
clusters observed with ROSAT (Eckert et al. 2012) 
at the outskirts (0.65-1.2r2oo)- 

• Within r5oo, the entropy radial profiles from the 
XMM-Newton and Suzaku observations follow a 
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Fig. 15. — Distributions of galaxies around Abell 1835 in a slice of 
photometric redshift, \z — z c \ < 0.0125. The FWHM for a Gaussian 
kernel is represented by the hatched circle at bottom right. The 
boxes represent the FoVs of the XIS pointings. The circle is the 
virial radius determined by lensing analysis. 
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Fig. 16. — Comparison of temperature deviations AT/Tah in 
the cluster outskirts (rsoo ~ r ~ r v i r ) and number density contrast 
8 = An/ (n) of large-scale structure for the region of (1— 2)r v ; r (left) 
and (2— 4)r v ; r (right). Circles (red) and diamonds (blue) represent 
the results for Abell 1835 and Abell 1689, respectively. Errors are 
the 68% CL uncertainty. Spearman's rank correlation coefficient, 
r s , gives that theses correlations do not occur by accident more 
than 90% (left) and 67% (right) confident. 

power-law form with a fixed index of 1.1 (K cx 
r 11 ), as predicted by models of accretion shock 
heating. In contrast, beyond r 50 o, the entropy pro- 
files become flatter in disagreement with the r 11 
relationship. 

• The H.E. and lensing mass estimates within rsoo, 
except within l.'O from the center, are consistent 
within errors (Figure 10). The H.E. mass pro- 
file unphysically decreases with radius at r £ r^oo . 
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TABLE 10 

Spearman's rank correlation coefficient between temperatures in 
cluster outskirts and density of large scale structure 



Name 


(1-2) 


^vir 


(2-4) 

r s a 


^vir 

pb 


Abcll 1835 
Abell 1689 
Two clusters 


0.800 ±0.122 
0-800+S 

0.762 ±0.133 


0.200 ±0.122 

200"*" ' 323 
u.zuu_ Q 200 

o.o28ig:°s 


0.600 ±0.171 
0.800t°;3°° 
0.571 ±0.177 


0.400 ±0.171 

200 +(1 ' 247 
u.zuu_ 200 

Q ,00+0.195 


a Spearman's rank correlation coefficient 
The probability of a null hypothesis 



Accordingly, the lensing masses are systematically 
higher than the H.E. masses in the cluster outskirts 
( ~ ^500 )• This means that most of the ICM in 
the cluster outskirts is out of hydrostatic equilib- 
rium, indicating additional pressure supports such 
as turbulence, bulk velocity and/or high ion tem- 
perature. 

• The gas mass fraction profile, combined with lens- 
ing and gas masses, agrees with ~ 90% of the 
cosmic mean baryon fraction from the WMAP 7- 
year results (Komatsu et al. 2011), in the range 
of 

^500 ~ r <-S r vh - In contrast, the H.E. based gas 
mass fraction profile, /g as (< r) = M gas /M H .E., 
continuously increases with radius beyond r.500, 
exceeding the cosmic mean value, as reported in 
Perseus cluster (Simionescu et al. 2011). These re- 
sults indicate that the breakdown of the strict hy- 
drostatic equilibrium, rather than the gas-clumping 
effect, is significant in the cluster outskirts. 

• The pressure profile inside the cluster agrees with 
that derived by stacked SZ flux for 62 nearby mas- 
sive clusters from the Planck survey, supports the 
reliability of Suzaku measurements of the thermal 
pressure. 

• The radial profiles of temperature, electron density, 
and entropy are independent of the halo concentra- 
tion of lensing NFW mass models for Abell 1835 
and Abell 1689, statistically corresponding cluster 
ages. It implies that the thermalization timescalc 
is much shorter than a difference of ages. The fur- 



ther statistical study with larger sample would be 
important to obtain more robust result. 

• The temperatures in the outskirts have azimuthal 
variation greater than measurement uncertainties, 
though the significant level is low. The electron 
density, and entropy profile in the South, West, 
and North directions are consistent within errors 
with one another. The temperature and entropy 
in East direction is higher than those in other re- 
gions, while the electron density is lower. We in- 
vestigate the correlation between the temperature 
in the outskirts (rsoo ^ r ^ r v i r ) and the large-scale 
structure, using Suzaku X-ray and SDSS photomet- 
ric data for Abell 1835 and Abell 1689. We found 
by Spearman's rank correlation coefficient a signif- 
icant correlation with the large-scale structure in 
(l-2)r v i r . The hot and cold temperature regions 
contact filamentary structure and low-density void 
regions outside the clusters, respectively. The cor- 
relation is not occurred by chance with more than 
90% confident. 
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man BMBF through the Verbundforschung under grant 
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APPENDIX 
POINT SOURCE ANALYSIS 

We would like to excise point sources because we are only interested in the ICM. As for the point-source subtraction, 
we first analyzed the XMM- Newton data (Observation ID=0 147330201) in which faint sources were resolved better 
than the Suzaku data. We detected 12 point sources using ewavelet task of the SAS software version 8.0.0 with a 
detection threshold set at 7er and used surround annular region for background subtraction. The source extraction 
radius is 30", and surrounding background ring in estimating the flux is defined by 30"-60", respectively. For the 
individual sources, we carried out spectral fits for the MOS1 and MOS2 spectra simultaneously using the same spectral 
model pegpowerlaw which offered photon index and flux in selected energy band. We fitted spectra in the energy range 
2.0-5.0 keV. We show the best-fit parameters for the individual point sources in Table 11. Those fluxes were higher 
than 2xl0 -14 erg cm~ 2 s _1 in the energy range 2-10 keV. 

We also searched for point sources with Suzaku, finding additional 20 sources outside the XMM-Newton FoV by 
CIAO tool wavdetect. We performed spectral fits to all the point sources with Suzaku according to the following 
procedure. As in the case of XMM-Newton, we jointly fitted the XISO, XIS1, and XIS3 spectra by pegpowerlaw. The 
source extraction radius is l.'O, and the NXB was subtracted before the fit. We fitted spectra in the energy range 
2.0-7.0 keV and excluded the point sources. The fluxes of these sources were higher than 2xl0 -14 erg cm -2 s _1 in 
the energy range 2-10 keV. 
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TABLE 11 

Best-fit parameters of detected point sources 



No. 1 



XMM-Newton(MOSl+MOS2) 
Photon Index Flux c Reduced-^ 2 (x 2 /d.o.f) 



Suzafai(XIS0+XISl+XIS3) 
Photon Index Flux c Reduced-x 2 (x 2 /d.o.f) 



1 


(210.365, 


2.935) 


— U.O 


40.0+^-jj 


0.74 (20.8/28) 


l-5 + 02 

— U.«I 


47.7+i?- 8 . 

— O.D 


0.78 (138.4/178) 


2 


(210.438, 


2.893) 




+ 25.7 

Zo.o_22 Q 


0.70 (16.2/23) 


y 2+0.3 

— 0.3 


27.1±t;I 


0.87 (92.6/106) 


3 


(210.323, 


2.731) 




tyr\ fj + 31.2 
— lo.O 


1.00 (23.9/24) 


— U.o 


14.8+H 

— d.Z 


n 7n 
U. f U 


f AQ £ /fiO 1 

(^o.D/oy ) 


4 


(210.221, 


3.018) 


1 3+0:9 

— 0.9 


20.7 + "„ 5 
— o.y 


0.57 (15.4/27) 


2.1+ 06 
— 0.6 


4.0+i-g 


i.Uo 


(Ad 1 /Qfl 1 ^ 

^4u.i/ oy j 


5 


(210.288, 


2.948) 


0.4+J-j 


16.2 +22 3 8 


1.07 (27.7/26) 


i.4+2-| 

— U.4 


14.5+^ 


n fifi 


^o4.4/ DZ J 


6 


(210.306, 


2.671) 


1.4 (fixed) 


14.9+|-0 
— o.u 


0.75 (12.8/17) 


— U.4 


8.5+11 


n fi/i 
U.D4 


(ZO.O/ 40 J 


7 


(210.291, 


2.727) 


1.4 (fixed) 


8.3+S-? 

— O.J. 


0.83 (20.0/24) 


l-5 + 08 
— U.o 


4.8+H 


1 1 fi 
1. lo 


(ol.o/ Z f J 


8 


(210.355, 


2.728) 


1.4 (fixed) 


6.7+6-* 


0.95 (19.9/21) 


l-8 + o 7 7 

— U. 1 




U. f 4 


/'OQ 7 /'SO' 1 ! 

i^ZO. / / oz ) 


9 


(210.140, 


2.971) 


1.4 (fixed) 


6.0+11 


0.82 (22.2/27) 


L1 +0.8 
— 0.8 


7-0+H 

— 0.0 


U.54 


/no a / r >Q\ 
(ZO.O/ ZO J 


10 


(210.251, 


2.955) 


1.4 (fixed) 


5.3+l : » 


1.15 (33.4/29) 


1 q+0.6 

-L.O (J g 


7 3+ 34 


n fin 


foo n /Qi^ 
^zz.y/ 00 J 


11 


(210.313, 


2.805) 


1.4 (fixed) 


4.3+i : ; 


1.32 (35.7/27) 


1.4 (fixed) 


7.9+o 1 

— Z.Q 


n o/i 
U.y4 


^ZO.O/ ZD J 


12 


(210.414, 


2.956) 


1.4 (fixed) 


1.8+^1 


0.87 (19.2/22) 


9 1+O.6 
5 


5.3+H 


n Qfi 
U.ou 


yoD. 0/ 41 ) 


13 


(210.424, 


3.016) 








o.9i°;« 


11.3+°+ 


n T7 


K Z{ . ( / 00 J 


11 


(210.554, 


2.832) 








1.4 (fixed) 


2.1+11 


l.zz 


^z4.4/ zU ) 


15 


(210.532, 


2.954) 










5-4+1-'? 


n fio 

U.OZ 


^14. z/ ZO ) 


16 


(210.495, 


2.988) 








l.o+i'i 


6.3+6+ 


n fin 


^14.0/ Z4 J 


17 


(210.544, 


2.864) 








i.6+°i 

— U.O 


8.1+fl 


n fin 
U.uy 


/lot l r i'7\ 

ylo.o/ 11 ) 


18 


(210.458, 


2.810) 








2.6+li 


3.1±|| 


n nfi 


fOA 1 /OC^ 

I ,£4.1/ zo J 


19 


(210.405, 


2.768) 








1.6+11 




0.66 


(13.2/20) 


20 


(210.380, 


2.754) 








— U.O 


W.6+li 

— o.z 


0.83 


(34.1/41) 


21 


(210.317, 


2.754) 








2 5+ - 7 
z -°-0.6 


c 7+2.0 
°-'-1.8 


0.98 


(38.4/39) 


22 


(210.296, 


2.561) 








2 6+ 10 


6.3+|-9 


0.95 


(22.7/24) 


23 


(209.969, 


2.902) 








2 4+ - 8 


3 4+ 2 ' 2 


0.69 


(24.3/35) 


21 


(210.184, 


2.771) 








1 8+ 1 - 1 
1 - 8 -0.9 


5 0+ 3 - 6 

O.U_ 2 g 


0.55 


(13.2/24) 


25 


(210.154, 


2.725) 








1 3+ 1 - 1 


r c+6.3 
J - d -3.6 


0.67 


(14.7/22) 


26 


(210.085, 


2.948) 








2 0+ 10 
z - u -0.9 


4 3+2-4 
*'°-2.0 


0.71 


(20.7/29) 


27 


(210.009, 


2.967) 








1.4 (fixed) 


6 5+ 3 - 7 
°-°-6.5 


0.74 


(15.5/21) 


28 


(210.049, 


2.751) 








2 1+ - 7 


6 2+ 2 ' 3 


1.16 


(44.2/38) 


29 


(210.161, 


3.129) 








1 6+ - 8 
i - D -0.8 


c 9+2-5 
°- z -2.1 


0.83 


(25.6/31) 


30 


(210.080, 


3.067) 








1.4 (fixed) 


11 Q+5-6 
ii - d -5.6 


0.80 


(16.8/21) 


31 


(210.083, 


3.094) 










r- 0+3.3 
O.o 2 7 


0.65 


(15.0/23) 


32 


(210.221, 


3.141) 








i-s±i:S 


2 5+ 1 ' 9 


1.07 


(23.6/22) 



a Serial number for point source. 

k Position of the point source. 

c The 2.0-10.0 keV flux in units of 10" 



erg cm 



In the observation area of Abell 1835 with Suzaku (~ 1000 arcmin 2 ), the number of point sources brighter than 
2xl0 -14 erg cm" 2 s _1 and l.OxlO -13 erg cm -2 s _1 expected from deep fields observations of XMM-Newton, Chandra, 
and ASCA (Kushino et al. 2002) are ^30 and ~3, respectively with T = 1.4. The actual numbers of point sources 
detected in this study above the two thresholds are 32 and 6 with T = 1.4, respectively, and are consistent with above 
expectations considering the Poisson error. 

We excised all the point sources detected in either the Suzaku or XMM-Newton observations. Normally we excluded 
a region of l.'O radius but used 2.'0 radius for a source (No.l in Table 11). Although we excluded circular regions 
around the point sources, the signals from brightest point sources are expected to escape from the excluded circular 
regions of 1 f radius, since the PSF of Suzaku is in HPD (Serlemitsos et al. 2007). Therefore, we simulated 
this residual signals from brightest point sources (No.l, 2, 3, and 5 in Table 11). Using the derived crude spectra of 
the 4 point sources, we simulated the XISO, XIS1, and XIS3 images using xissim for the observation that each point 
source is present. We then subtracted the residual signals for all annular spectra before the spectral fitting. In order 
to prevent the deterioration of the statistics, we only considered the residual signals from brightest four point sources. 
The Abell 1835 field of view contains two low- mass groups (see e.g. Bonamente et al. 2012). One was excised in our 
analysis (No. 21 in Table 11). Another was not because its flux was less than 5% of the ICM flux in the corresponding 
annulus (6.'0-9.'0). Consequently, the spectral analysis did not suffer from the emission from these two groups. 

6'-9' AND 9'-12' REGIONS IN THE EAST DIRECTION 

It was difficult to constrain the model parameters in the 9'0— 12f0 region of the East direction because of the influence 
of bright point sources (see Figure 1). We could not constrain the 90% confidence level upper limit of the temperature 
nor calculate deprojected electron density for the 9'0— 12^0 region in the East direction (see Table 6 and 7). Therefore, 
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TABLE 12 

Best-fit values of the 6.'0-9.'0 and 9.'0-12f0 regions 

FITTING IN THE EAST DIRECTION 



Case 


kT 
(keV) 


Norm a 


Reduced-^ 2 b 
(x 2 /d.o.f) 


GAL 

GAL+ICM 
GAL2 


a.^i_ 118 
S 49+2.25 

■ 4.76tt;| 


, 7t -+0.54 

2 16+ - 52 
1 7 4+0.49 
1 - '^-0.45 


1.26 (258/204) 
1.24 (253/204) 
1.30 (266/204) 



a Normalization of the apec component scaled with a factor of 
SOURCE.RATIO.REG/n e 

Norm = SOURCE.RATIO.REG J / [i7T (1 + z f D ^ ] x 

10 20 cm _J arcmin -2 , where Da is the angular distance to the 
source. 

k X 2 °f the fit when parameters are tied between the 6.'0— 9.'0 and 
9.'0-12.'0 regions in the East direction. 

Case-GAL9'-12'East „ Case-GAL+ICM 9'-12* East nl Case-GAL2 9'-12' Easi 




Energy (keV) Energy (keV) Energy (keV) 



Fig. 17. — The NXB-subtracted spectra of XIS3 (black crosses) and XIS1 (red crosses) for the 9.'0-12f0 region in the East direction fitted 
with the ICM model plus the X-ray background model described in Section 3.1. Left, center, and right panels correspond to the Case-GAL, 
Case-GAL+ICM, and Case-GAL2, respectively. The ICM, CXB, LHB, MWH1, and MWH2 emissions for the XIS3 spectra are shown in 
magenta, blue, orange, green, and cyan lines, respectively. Sum of the CXB, LHB, and MWH1 emissions for the XIS3 spectra are indicated 
by green-cyan line. The total model spectra of XIS3 and XIS1 are shown in black and red lines, respectively. The lower panels show the 
residuals in units of a. 

we fitted the spectra for the 6.'0-9'0 and 9f0-12.'0 regions in the East direction linked the ICM temperature and 
normalization in these regions. Figure 17 show the results of the spectral fit in the 9.'0-12f0 region for the three 
background cases. The best-fit parameters of ICM components and x 2 values in this region are listed in Table 12. The 
best-fit parameters were consistent within the systematic errors for the two regions. 
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